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The ground state of odd-A^ polyynic oligomers CjvH2 features kink-solitons in carbon-carbon 
bond-length alternation (BLA) patterns. We perform a systematic first-principles computational 
study of neutral and singly-charged kinks in long oligomers addressing relationships between BLA 
patterns, electron energy gaps, and accompanying distributions of spin and charge densities, both 
in vacuum and in the screening solvent environment. A quantitative comparison is made of the 
results derived with four different ab initio methods: from pure DFT to pure Hartree-Fock (IfF) 
and including two popular hybrid density functionals, B3LYP and BHandHLYP. A clear correlation 
is demonstrated between the derived spatial extent of kinks and the amount of HF exchange used 
in the functionals. For charged kinks, we find a substantial difference in the behavior of charge and 
spin densities. 
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I. INTRODUCTION 

One-dimensional (ID) semiconductor nanostructures 
such as TT— conjugated polymers (CPs), nanotubes and 
nanowires attract a great deal of attention being both 
interesting scientifically and important for technological 
advances. The nature and properties of excess charge 
carriers and electronic excitations on these structures are 
fundamental for many applications in (opto)electronics, 
energy harvesting and sensors. As covalent bonds in 
these systems provide for wide electronic bands, excess 
charge carriers have a high propensity for delocalization 
and a potential for high mobilities. The interaction of 
electrons with other subsystems such as displacements of 
the underlying atomic lattice and polarization of the sur- 
rounding medium can however put inherent limitations 
on the mobility and may even lead to self-localization of 
excess carriers. The electron-lattice coupling has been a 
subject of a great body of studies for various ID systems, 
and in particular for CPs (see, e.g., Refs. [U, Sjlll for re- 
views and original references). It is well-known that, as 
a result of the interaction with lattice distortions, elec- 
trons and holes in CPs may undergo self-trapping into 
polarons accompanied by localized bond-length modula- 
tion patterns and new features in the optical absorption 
due to local intragap electronic levels. Kink-solitons are 
another and qualitatively different type of localized in- 
trinsic "defects", which are specific to CPs with the de- 
generate ground state and exhibit topological kink-like 
patterns in bond-length modulations (Fig. [T|) and even 
deeper lying intragap levels. Solitons are believed to play 
a prominent role in the physics of trans-polyacetylene (for 
a recent discussion, see Ref. Q). Polaronic states, being 
non-topological excitations, are generic for a broad range 
of CPs. Polarons in CPs can also be thought of as bound 
kink-anti-kink pairs. 

Numerous first-principles calculations at various the- 
oretical levels have been applied to study the struc- 
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FIG. 1: Carbon-carbon bond length alternation patterns for 
neutral triplet kink-solitons on polyynic oligomers in vacuum. 
The circular data points show results derived from computa- 
tions with end effects eliminated, solid lines depict analytical 
fits as discussed in the text. The computational methods used 
are indicated next to the respective curves. 



ture of the ground state and excitations in CPs. With 
such calculations one expects to be able to elucidate ef- 
fects due to valence electrons and many-electron inter- 
actions. Interestingly enough, despite a relatively long 
history of electron-lattice polarons in CPs (see, e.g., mul- 
tiple references in Refs. [1,0])) certain questions have 
been raised recently regarding applicability of different 
ah initio frameworks to describe the polaron formation. 
That pertains to the reported failure of the local-density- 
approximation and generalized-gradient-approximation 
DFT (density-functional theory) schemes to detect self- 
localized charge density distributions in various charged 
oligomers, while Hartee-Fock (HF), parameterized semi- 
empirical and possibly hybrid-functional DFT ,7| meth- 
ods have been reported to lead to charge localization in 
the middle of oligomers (e.g., Refs. i,i,i,[i3). Quite 
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different level-of-theory-dependent magnitudes of the ef- 
fective electron-phonon coupling have also been found for 
the dimerized ground states of polymeric systems [lTI.[l^. 
These observations thus bear on an important general is- 
sue of a choice of appropriate ab initio methods 0, 
to faithfully describe properties of ID semiconductors. 

In this paper we contribute to this discussion by com- 
paring results on the structure of kink-solitons in long 
odd- TV polyynic oligomers C ArH2 obtained with different 
ab initio methods including pure HF and DFT (PBE) 
schemes as well as two popular hybrid-functional frame- 
works (B3LYP and BHandHLYP). Structurally simple 
polyynic linear carbon chains naturally lend themselves 
to be a prototypical example of ID semiconductors. It 
should be noted however that polyynic chains continue 
to be the subject of much attention in their own right 
[H [lljlljli, m, [13 • Interestingly, they were pre- 
dicted [18l.ll9j| to possess a rich family of electron-lattice 
self-localized excitations because of the extra spatial de- 
generacy of molecular orbital levels. Neutral kinks take a 
special, in a sense, place representing in fact the ground- 
state structure of odd-A^ oligomers. By studying "long 
enough" oligomers, we are trying to find the "natural" 
spatial extent of the kinks as determined by inherent sys- 
tem interactions rather than by the end effects due to the 
hydrogen termination. As Fig. [1] illustrates, the spatial 
extent of solitons turns out to be very strongly depen- 
dent on the computational method used and evidently 
correlated with the magnitude of the equilibrium lattice 
dimerization. 

An addition of an electron (or a hole) to oligomers re- 
sults in charged kinks and increases their spatial extent. 
The strong dependence on the computational method 
used is also reflected in spatial distributions of charge and 
spin densities associated with the kinks. Along with the 
effects of electron-lattice interactions, in this paper we 
also study how the surrounding polar solvent influences 
the resulting structures, a much less explored topic. The 
solvent environments are ubiquitous in many applications 
involving fundamental redox processes. While neutral 
oligomers are practically not affected by the surrounding 
polar medium, the charged kinks are found to strongly 
interact with the solvent. We have recently demonstrated 
[20I [2l| at the level of ab initio computations that solva- 
tion of excess charge carriers on even-A^ oligomers results 
in the formation of polaronic states, and that solvation 
and lattice distortions can act synergistically to increase 
the degree of charge localization. As the kink structure 
is present already on neutral odd-A'^ oligomers, one can- 
not characterize it as resulting from self-localization of 
an excess electron. We however observe that the spa- 
tial extent of the charged kink lattice distortion pattern 
gets shorter upon solvation, and in this sense the syner- 
getic effect also takes place. This qualitative observation 
is independent of the computational method used, as is 
the fact that the solvent environment affects the spatial 
distribution of charge densities much more significantly 
than that of spin densities. This may be an indication of 



a different character of many-electron system responses 
to charge and spin perturbations. 



II. COMPUTATIONAL METHODS AND DATA 
PROCESSING 

All ab initio computations in this study were per- 
formed using the GAUSSIAN 03 suite of programs [23 |. 
We have employed purely Hartree-Fock (HF) and DFT 
(PBE) calculations along with the well-known hybrid- 
DFT functional BHandHLYP and B3LYP in order to 
compare results of the different levels of theory. We re- 
mind the reader that the amounts of HF exchange in- 
cluded in these hybrid methods are 50% for BHandHLYP 
and 20% for B3LYP. A rich 6-3U++G{d,p) all electron 
basis set has been employed throughout. It should be 
noted that wherever possible, the results of our "in vac- 
uum" calculations have been verified to compare well 
against previously published ab initio data on even-A^ 
[H [III and odd-A^ [H ^ neutral and charged CatHs 
systems both in terms of energetics and optimized bond 
lengths. In order to investigate the effects of solvation 
in a polar medium, GAUSSIAN 03 ^22^ offers its imple- 
mentation of the polarizable continuum model (PCM) 
described in original publications [H, [H, [H, . For all 
"in solvent" PCM calculations, water has been chosen as 
a polar solvent with its default parameters in GAUSSIAN 
03. 

In working with long oligomers we met with certain 
computational method-dependent limitations on achiev- 
ing satisfactory levels of convergence and spin contami- 
nation. In order to produce reliable results, therefore, we 
had to resort to using both restricted (RO, for HF and 
BHandHLYPJand unrestricted (U, for B3LYP and PBE) 
calculations [2^. The longest odd-A^ oligomers studied 
have been A^ = 41 (and some A^ = 61) for HF, A^ = 61 
for BHandHLYP, N = 101 for B3LYP, and A^ = 141 for 
PBE methods. To our knowledge, these oligomer lengths 
are appreciably longer than used previously in compara- 
ble calculations. 

The physical quantities of interest to us in this pa- 
per are the equilibrium (relaxed) geometries of relative 
carbon atom positions, atomic charge densities (to be 
denoted Qn) and atomic spin densities (pj). We refer 
the reader to our previous publication [2l| for an exten- 
sive discussion of some aspects of a retrieval of charge 
densities. The charge densities g„ displayed below in 
Fig. [5] are excess charges understood as the difference of 
the raw charges in charged- and neutral-kink configura- 
tions. The raw atomic charges, in turn, are derived from 
Lowdin population analysis [12, [1^. We note that, for 
excess charge, the difference between Lowdin and Mul- 
liken [13, |2§| populations is relatively small ^l| . For the 
atomic spin densities p^, shown in Figs. [S] and [71 we use 
the only available MuUiken population data. 

A very convenient and well-known P, |3, [ll|, [l2| way to 
characterize the geometry of dimerized polymeric struc- 
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FIG. 2: BLA patterns for a series of increasing chain length 
neutral odd-TV oligomers (from A*' = 21 to = 101) com- 
puted using the B3LYP method, empty data points. Filled 
data points show results for the C21H2 oligomer derived in 
Ref. 

tures is via bond-length alternation (BLA) patterns 

^ {-ir {In - In-l) , (1) 

where In is the length of the nth carbon- carbon bond 
defined, e. g., as on the right of the nth carbon atom 
(we do not discuss the end bonds to hydrogens). For 
representation in this paper, we chose the index of the 
central carbon atom of the oligomer as n = 0. In the in- 
finite dimerized neutral structure, the dimerization pat- 
tern ([J) would be uniform, that is, independent of the 
spatial bond position n, we denote the magnitude of this 
equilibrium ground-state pattern as Sq. This is not to 
be confused with a specific site value of Sn=o, the lat- 
ter is not being explicitly used anywhere in the paper. 
The double degeneracy of the ground state of the infinite 
polymer is in that the uniform pattern ([1]) can assume 
either value of +So or —So [Ull]. For our purposes, we re- 
trieve the magnitude from the middle segment of long 
even-A^ oligomers featuring a well-established spatially- 
independent behavior. 

As an example. Fig. [2] displays the relaxed non- uniform 
BLA patterns for a series of neutral oligomers CArH2 
(A^ = 21, 41, 61, 81, and 101), all obtained within 
B3LYP computations. Along with our results also shown 
and in very good agreement is the data from Ref. [l3| 
on C21H2 oligomer (the longest they studied was with 
N = 23). The figure makes it transparent that the kink 
structure continues to evolve in longer oligomers well be- 
yond A^ = 21, and the latter case cannot serve even ap- 
proximately to estimate the inherent extent of the kink. 
While illustrating the importance of a sufficient length 
of the oligomer, the figure also clearly shows that the 
end effects on the BLA pattern are quite substantial and 
persistent. The "shape" of the end-specific BLA struc- 
tures stabilizes and becomes easily discernible in longer 
oligomers, it, of course, coincides with the one observed 
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FIG. 3: Illustrating end effects in BLA patterns and their 
elimination: Panel (a) compares the neutral kink structure 
on C101II2 (solid line) with the BLA patterns of "halves" of 
the neutral C100H2 oligomer (dashed), both obtained with 
the B3LYP method. After the subtraction of the common 
end BLA structures in this case, the kink BLA pattern of 
panel (b) results. Panel (c) compares relaxed BLA patterns on 
neutral C101H2 and C141H2 (data points) as well as "halves" 
of BLA patterns of the ground state of C140H2 (dashed lines), 
all obtained with the PBE method. 

in BLA patterns of the ground state of even-A^ oligomers. 
This fact allows to "eliminate" the end effects by com- 
parison and an appropriate "subtraction" procedure, as 
illustrated in Fig. [3{a) and (b) . Such a procedure of re- 
moving the end effects from the finite oligomer data has 
been used to produced the kink BLA patterns shown in 
Fig. [1] and other figures in this paper. 

III. NEUTRAL KINKS 

To quantify the spatial extent of kink-solitons dis- 
played in in Fig. [1] one can conveniently use the ana- 
lytical shape of 

n 

Sn = Sq tanh y (2) 

well-known [l|, from the continuum electron-phonon 
models of CPs. Since the ground-state BLA magnitude 
Sq has already been determined from the computation 
of long even-A^ oligomers, the soliton half-width I is the 
only parameter to be found from the fitting to the ac- 
tual data. As n in Eq. ([2|) is the dimensionless index, 
parameter I measures the spatial extent in units of the 
corresponding carbon-carbon lattice constants. The re- 
sulting fits for each of the computational methods used 
are shown in Fig.[T]and appear quite representative. The 
numerical values of the optimal parameters are collected 
in Table [J The table also features the magnitude Eg of 
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TABLE I: Numerical values obtained for the dimerization 
magnitude Sq, neutral kink half- width I, and HOMO-LUMO 
gap Eg with different methods from the longest oligomers 
studied. 



method 


So (A) 


I 


Eg (eV) 


ROHF 


0.183 


1.5 


8.59 


ROBHandHLYP 


0.135 


3.5 


4.18 


UB3LYP 


0.091 


10.8 


1.77 


UPBE 


0.036 


fa 29 


0.41 
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FIG. 4: The dimerization magnitude (So and HOMO-LUMO 
gap Eg as functions of the inverse neutral kink half- width 1 // 
obtained with different methods, see Table |T] Lines connect- 
ing data points are used here only to guide the eye. 



the gap between the highest occupied (HOMO) and low- 
est unoccupied (LUMO) molecular orbitals of the neu- 
tral even-A^ oligomers. The values presented in the table 
have been derived from the data on the longest oligomers 
studied: I for the kink structures in Fig. [l] 5o and Eg for 
N = 100 oligomers in the case of HF, HandHLYP and 
B3LYP methods, and for N = 140 in the case of PBE. 

We should note here that a careful examination of the 
fitting to the PBE-derived kink BLA pattern on C141H2 
suggests that even on such a long oligomer, there is still 
some "pressure" exerted on the kink by the oligomer 
ends, and the derived value of / ~ 29 may be a bit under- 
estimated. At the same time, there are clear indications 
that a fully-relaxed kink structure with a self-consistent 
spatial size can be found formed within a purely DFT- 
type computational framework. We consider this a useful 
observation bearing on the inherent qualitative capabil- 
ities present in such frameworks, even if the numerical 
magnitudes may be off the mark. In the other extreme of 
the HF computations, the resul ting kink width (in agree- 
ment with earlier calculations [23|) is so small that the 
continuum-like analytical form ^ also has its limitations 
in representing some of the data points. 

The overall picture of the relationships between results 
obtained with different methods is displayed in Fig. 3] in 



the form of Sq and Eg as functions of the inverse kink 
half-width l/l. This view is motivated by the continuum 
electron-phonon models of CPs ([H, [1] and, particularly 
for polyynes, where the kink energy gap param- 

eter 



A(a;) ~ Ao tanh 



^0 



(3) 



as a function of the coordinate x along the (infinite) poly- 
mer depends on the equilibrium magnitude Ag and cor- 
relation length ^0 that are related as 



^0 = wf/Aq 



(4) 



(vp being the Fermi velocity). More generally, relation- 
ships like Eqs. ([3]) and ^ are also found for kink-solitons 
of nonlinear field theories, such as cj)'^ [30] . The HOMO- 
LUMO gap Eg in such models is just 2Ao, and our di- 
mensionless half- width / oc ^0, hence the idea to see if the 
relationship 



Eg cx l/l 



(5) 



might work for our data. In purely electron-phonon mod- 
els (no electron-electron interactions), however, the mag- 
nitude of the gap parameter and the dimerization ampli- 
tude are just proportional to each other: Aq cx Sq, and 
one therefore could as well be verifying the relationship 



6q oc l/l. 



(6) 



Both the dimerization amplitude 60 and the kink spa- 
tial extent I reflect the strength of the effective electron- 
phonon interaction: the stronger the interaction, the 
larger 60 and the smaller I are expected to be. The 
data presented in Table |T] and in Fig. |4] clearly and con- 
sistently demonstrate that the effective electron-phonon 
interaction as revealed in results of ab initio computa- 
tions becomes stronger with the increased amount of HF 
exchange present in the energy functionals of the corre- 
sponding frameworks. Of course, the HOMO-LUMO gap 
Eg also shows a steep growth with the amount of HF ex- 
change. One should however recall that electron-electron 
interactions affect Eg not only "indirectly" via the effec- 
tive magnitude of the electron-phonon interaction (and 
resulting dimerization Sq) but also "directly" (as is easy to 
verify by keeping Sq artificially fixed at non-optimal val- 
ues). It is in this connection that it becomes interesting 
to compare our calculated data with both relationships 
and ([S]), that would be equally valid in non-interacting 
electron models. 

Already a cursory look at Fig. |4] is sufficient to see 
that deviations from simple proportionality ^ are quite 
substantial. On the other hand, the data is relatively 
much closer to relationship ([S]), although by no means 
exactly following it. We cannot exclude that this in- 
teresting observation bears on the possibility of (rough) 
mapping between results of different models. One inter- 
pretation of the observation is that it is indeed the "full" 
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FIG. 5: The spatial distribution of BLA patterns and the 
spin density in neutral triplet kinks as derived with difTer- 
ent methods: (a) ROHF, (b) ROHandHLYP, (c) UB3LYR 
The results are shown in forms scaled to unit maximum val- 
ues. The scaled raw (oscillating) atomic spin densities are 
displayed as simple solid lines, their averaged (as described in 
the text) behavior is indicated with data symbols. 



energy gap parameter ^ (rather than just the lattice 
displacements) that is a (more) adequate order parame- 
ter to describe spatially inhomogeneous structures within 
an effective nonlinear field model. 

The polyynic structures, in comparison with poly- 
acetylene, possess an extra degeneracy of 7r-electron 
single-electron levels, resulting from two independent ori- 
entations of wave functions in the plane perpendicular to 
the polymer axis. The midgap electronic level formed in 
the presence of the kink-shaped BLA pattern has there- 
fore an internal degeneracy of four (including spin) and 
is occupied by two electrons when the kink is uncharged 



0, [T^ . The neutral kink therefore can be either in the 
singlet or triplet states. In the non-interacting electron 
models, these states are degenerate. It is known though 
from the previous studies [13, that electron interac- 
tions prefer the triplet state. In full agreement, all our 
simulations have also resulted in the triplet, spin S" = 1, 
ground state of neutral kinks. Only these triplet states 
are being explicitly presented in this paper. In the con- 
text of our discussion, this comes in "handy" as one can 
compare the spatial structures of BLA patterns with non- 
vanishing spin density distributions. 

To facilitate this comparison, in Fig. [5] we plot BLA 
patterns and spin densities (retrieved from the MuUiken 
population analysis) in the scaled form emphasizing their 
spatial structures. (From now on, comparison results 
would be shown as derived with three methods exclud- 
ing PBE.) The raw atomic densities exhibit character- 
istic oscillations from carbon to the next carbon. To 
better see the spatial extent of the spin density. Fig. [5] 
also displays the averaged behavior of scaled atomic den- 
sities, where a simple-minded nearest-neighbor averag- 
ing of site-n specific quantities r„ is done as (r„) = 
0.5(r„ -I- 0.5(r„_i + r„+i)). (We use the same type of 
averaging for data in Figs. [6] and [7l) It is well known [2^ 
that calculations of atomic charges from ah initio results 
depend on the basis set transformations and sometimes 
lead to artifacts and spurious results. While we cannot 
exclude that some artifacts may be reflected in the raw 
atomic-centric data, we expect the averaged behavior to 
be much less susceptible to such dependencies (compare 
to our discussion [21| of atomic- and cell-centric charge 
densities). With the averaged data, it is especially trans- 
parent from Fig. [5] that for all computational methods 
used, there is a very good correspondence of the spatial 
extent of the kink BLA structure and the region of the 
spin localization. Needless to emphasize that, of course, 
each of the methods yields a very different spatial size 
of this region. In terms of the magnitude of oscillations 
of atomic densities, the panel (c) of course stands out, 
this may be related to the fact that the data for that 
panel was generated with an unrestricted computational 
method, while panels (a) and (b) have been produced 
with restricted computations. 

As was expected from our earlier analysis [13, HH of 
neutral even-A'^ oligomers, we have found practically no 
effect produced by solvation on the properties of neutral 
kinks at all levels of computation. That is why the data 
shown for neutral kinks has been restricted to results "in 
vacuum" only. The situation changes substantially for 
charged kinks. 



IV. CHARGED KINKS AND THE EFFECTS OF 
SOLVATION 

In this Section we discuss characteristics of singly- 
charged kinks resulting upon the addition of an extra 
electron to long odd-A^ polyynic oligomers. (These sys- 
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FIG. 6: BLA patterns and excess charge distributions for singly-charged kinks as obtained with different methods: ROHF in 
columns (a), ROHandHLYP in (b), and UB3LYP in (c). The first row displays results for chains in vacuum and the second for 
chains in the solvent environment. To compare, BLA structures for neutral kinks are also depicted with dotted lines. The raw 
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FIG. 7; As in Fig. |6] but for spin (instead of charge) density distributions. 



terns and nearly charge-conjugation symmetric and sim- 
ilar results are obtained upon the addition of an extra 
hole [21].) Our discussion will revolve around the data 
shown in Figs. [6] and [7] displaying resulting equilibrium 
distributions of BLA patterns, excess charge and spin 
densities. The overall charge of the system is of course 



just the electronic charge and the spin S — 1/2. While 
different computational methods yield numerically differ- 
ent magnitudes, common qualitative features are clearly 
discernible in the figures. 

All results indicate localized characteristics around the 
center of the kink in the middle of oligomers. This should 
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TABLE II: Electron affinities in eV calculated with various 
methods in vacuum and in the solvent for singly-charged 
oligomers of different lengths. 



method 


C41 


Ha 




H2 


Cioi 


H2 




vac 


solv 


vac 


solv 


vac 


solv 


ROHF 


2.39 


3.66 


2.43 


3.66 






ROBHandHLYP 


3.69 


4.68 


3.78 


4.69 






UB3LYP 


3.88 


4.68 


4.09 


4.72 


4.23 


4.73 



not be construed as self-localization of an excess electron 
since midgap levels are formed already within the neutral 
kink structure and can be readily occupied by that elec- 
tron. As perceived via BLA patterns, charged kinks evi- 
dently become wider in their spatial extent than neutral 
kinks, which would not be happening in non-interacting 
electron models [IS, JLSf] . This appears natural to ratio- 
nalize as the effect of Coulomb forces trying to spread 
out the excess charge, as if there was an effective width- 
dependent capacitance associated with the kink. This 
view is further supported by the observation of an ap- 
preciable effect of the screening solvent environment on 
the spatial extent of the kink. As is clearly seen, solvation 
decreases the extent of charged kinks in all computations 
and practically restores the neutral kink BLA pattern in 
the case of B3LYP (columns (c)). 

We did not attempt to analyze the length of charged 
kinks by fitting to the analytical dependence ([2]) as even 
longer oligomers would need to be studied to avoid the 
influence of end effects. In addition to the visual per- 
ception of those in the figures, the data in Table |TT] illus- 
trates oligomer-length-dependent progression of the elec- 
tron affinities EA understood as the difference of the total 
(relaxed) system energies of neutral and singly-charged 
kinks: 



One can see that this quantity in the vacuum data has 
not converged yet to self-consistent (length-independent) 
values. In contrast, the convergence has been (nearly) 
achieved for the "in solvent" data. That corresponds 
very well with the picture of a great effect of the solvent 
on the spread of the excess charge density in Fig. [H] which 
is obviously much farther away from oligomer ends in the 
solvated case (second row). Comparing EAs in vacuum 
and in solvent, the table also shows the energetic signifi- 
cance of solvation: from more than 1.2 eV in HE compu- 
tations to 0.5 eV in B3LYP. The method dependence here 
is consistent with the fact that the spatial extent of the 
excess charge is shortest in HE, hence the electrostatic 
Coulomb effects are strongest. 

A remarkable and consistent difference is transparent 
in the behavior of charge and spin densities in Eigs. [5] 



and [7] for all computational methods. The spatial spread 
of the charge densities in vacuum is much more signif- 
icant than of the spin densities. The former strongly 
respond to the solvent environment, while there is very 
little effect on the latter. Even in the solvent the charge 
and spin densities appear perceptibly distinct. This dif- 
ference is at variance with what one would expect in a 
simple picture of a single localized electron, whose prob- 
ability density would establish both charge and spin den- 
sities. Evidently, it is the response of the many-electron 
system that should account for the observed difference. 

Given the fact that the solvent so strongly affects the 
excess charge distribution, one could relate the effect to 
long-range electrostatic Coulomb forces that polarize the 
many-electron system of the polyyne and cause remote 
variations of the charge density but not the spin density. 
This however does not appear as a complete explanation 
as the difference is observed for the averaged densities as 
well ("true" excess charge). Apparently, there may be 
inherent distinctions in system responses here to charge 
and spin perturbations, perhaps related to the nature of 
the order parameter. While, of course, not directly appli- 
cable to our system, one cannot help speculative wonder- 
ing if there could be some relationship to the phenomenon 
of spin-charge separation well-known in ID electronic liq- 
uids [m. 

At the same time, we recognize that some computa- 
tional artifacts may have affected the observed picture. 
We note that the overall spatial behavior of spin densi- 
ties on charged kinks is similar to observed on neutral 
kinks, Eig.O and one is reminded that it is only B3LYP 
computations (column (c) in Eigs. [6] and [7]) that have 
been performed in the unrestricted mode. This may have 
contributed to larger degrees of spatial variations of the 
atomic spin densities observed in that case. It is not clear 
if more significant spin responses would have resulted 
from unrestricted calculations with HE and HandHLYP 
functionals. 



V. DISCUSSION 

Structural simplicity of polyynic carbon chains CArIl2 
makes them a good prototypical example of ID semicon- 
ductors and a convenient "playground" for using vari- 
ous theoretical methods. In this paper we extended the 
previous work on these systems to systematically study 
characteristics of kink-solitons that are naturally formed 
in bond-length alternation patterns of odd- TV oligomers. 
We attempted a quantitative comparison of results that 
are obtained at different well-known levels of ab ini- 
tio computations: from pure DET (PBE) via hybrid 
B3LYP and BHandHLYP to pure Hartree-Eock (HE). 
These methods have been listed here in the order of the 
increased amount of HE exchange in their functionals (0, 
20, 50, and 100%, respectively). The earlier findings of 
the increased magnitude of the uniform lattice dimer- 
ization pattern (Sq) and the HOMO-LUMO gap (Eg) 
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with the increased amount of HF have now been comple- 
mented by expHcit observations of the decreased spatial 
extent (l) of the kinks. This is consistent with stronger 
effective electron-phonon couplings "induced" by HF ex- 
change. By studying long oligomers, we tried to min- 
imize the influence of the end effects and to establish 
self-consistent sizes I of neutral kinks as determined by 
inherent system interactions. While numerical values of 
resulting parameters vary substantially from method to 
method, the data suggests that the relationship Eg ccl/l 
could serve as a rough first approximation for mapping 
results between different methods. This may be indica- 
tive of the energy gap parameter (including effects from 
both electron-phonon and electron-electron interactions) 
being an adequate order parameter for building a unified 
effective nonlinear field model. 

The more important issue, however, is if any (single) of 
these methods is adequate to describe real systems, which 
cannot be resolved without verification against experi- 
ments. Based on comparisons with structural an d op tical 
experimental data, it was suggested previously [ll| that 
BHandHLYP may be the "right" method to compute the 
equilibrium lattice configurations of uniform systems but 
B3LYP should be used for the electronic structure. The 
latter conclusion was motivated by BHandHLYP-derived 
HOMO-LUMO gaps Eg being larger than experimental 
"optical gaps" and would be unfortunate in the sense of 
the internal consistency of a single theoretical method. In 
view of the strong excitonic effects in fD systems [islls^. 
however, the situation may need to be reassessed to in- 
clude those effects in calculations of the optical spectra. 
On the other hand, the problem of small DFT-derived en- 
ergy gaps Eg has been successfully remedied [l^l by cal- 



culating corrections to the quasi-particle energies within 
the framework of the GW-approximation 0, [Tj, Hfi] (ex- 
citons are further calculated via the Bethe-Salpeter equa- 
tion \\S\). Even from a purely theoretical viewpoint, it 
would be quite interesting to compare results derivable 
for polyyncs with the help of GW to the ones obtained 
with less computationally intensive hybrid-DFT methods 
used in this paper. The picture emerged in our discus- 
sion of various computational outputs (Fig. [4]) strongly 
suggests, however, that they are intimately related, and 
just taking into account GW-corrections for Eg may turn 
out not to be sufficient, it is the optimal lattice structure 
of the ground state itself that could undergo changes if 
treated on the same footing. Where in Fig. [?] would then 
the corresponding data points fall appears a challenging 
question to answer. 

Consistently with our previous studies [l^, [2lj of po- 
laronic states and qualitatively method- independent, we 
have found a substantial localizing effect that the sur- 
rounding solvent produces on the excess charge densities 
associated with charged kinks. As a result of solvation, 
the equilibrium spatial extent of charged kinks also be- 
comes shorter; one may say that the lattice deformation 
and polarization of the solvent act synergistically to in- 
crease the degree of charge localization. We expect the 
drag from the solvent to significantly reduce the mobil- 
ity of charged kinks as it does for solvation-induced po- 
larons [ssl. |34|. Quite different responses in charge and 
spin densities found in all computations clearly indicate 
that electron-electron interactions are strongly involved 
in accommodation of an excess charge carrier by the kink- 
soliton structure. 
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